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Abstract. A focus of recent experimental and theoretical studies on heavy fermion systems close to antiferromagnetic (AFM) 
quantum critical points (QCP) is directed toward revealing the nature of the fixed point, i.e., whether it is an itinerant 
antiferromagnet [spin density wave (SDW)] type or a locally-critical fixed point. The relevance of the local QCP was proposed 
to explain the £/r-scaling with an anomalous exponent observed for the AFM QCP of CeCus qAuq.i. In this work, we have 
investigated an AFM QCP of another archetypal heavy fermion system Ce(Rui_;(;Rh;c)2Si2 with x = and 0.03 (~ x^) using 
single-crystalline neutron scattering. Accurate measurements of the dynamical susceptibility lmx{Q^E) at the AFM wave 
vector Q = 0.35c* have shown that lmx{Q^E) is well described by a Lorentzian and its energy width T{Q), i.e., the inverse 
correlation time depends on temperature as T[Q) = ci +C2r^/^^^-^, where ci and C2 are x dependent constants, in low 
temperature ranges. This critical exponent 3/2 proves that the QCP is controlled by the SDW QCP in three space dimensions 
studied by the renormalization group and self-consistent renormalization theories. 



INTRODUCTION 

Quantum critical points (QCP) separating ferromag- 
netic or antiferromagnetic states from paramagnetic 
Fermi liquid states in strongly correlated electron sys- 
tems have been investigated for decades. Successful de- 
scriptions of quantum critical behavior were provided by 
the self-consistent renormalization (SCR) theory of spin 
fluctuations |1, 2] for J-electron systems based on the 
Hubbard model. The mean-field type approximations in 
this theon^ were justified by the renormalization group 
studies ISlHl using Hertz's effective action above upper 
critical dimensions. For the ferromagnetic QCP, theoret- 
ical predictions were supported by experimental studies 
of J-electron systems 1 1 ] . In contrast there is little ex- 
erimental understanding of the antiferromagnetic QCP 
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A recent intriguing issue of QCP under controversial 
debate is directed toward revealing relevant fixed points 
for antiferromagnetic QCPs in heavy-fermion systems 
For energy scales much lower than the Kondo tem- 
perature Tk, / and conduction electrons form composite 
quasiparticles with a large mass renormalization in para- 
magnetic heavy fermions. By tuning a certain parame- 
ter, e.g., pressure or concentration, an antiferromagnetic 



long range order emerges from the Fermi liquid state. In 
a weak coupling picture, it has been hypothesized that 
the same QCP as the J-electron itinerant antiferromag- 
net, referred to as spin density wave (SDW) type QCP, is 
relevant to the heavy fermion QCP |6, 5]. 

However despite a number of experimental studies of 
heavy-fermion systems showing non-Fermi liquid be- 
havior, none of them definitively supports this QCP 
1 7, 8, 9, 10]. This stems partly from experimental dif- 
ficulty in measuring weakly divergent quantities around 
a QCP, especially for bulk properties, which has been 
also the case for J-electron itinerant antiferromagnets 
1 2]. On the contrary, several recent experiments suggest 
the possibility of a novel strong coupling picture of the 
QCP jslliillij]. Among these studies, single-crystalline 
neutron scattering investigations of the heavy fermion 
CeCu5.9Auo.i provided interesting insight 1111 . On the 
basis of the observed E / J-scaling with an anomalous ex- 
ponent 111], and effective two space dimensions flS*], a 
scenario of a locally critical QCP was proposed |5, 14]. 

In this work, we have studied the antiferro- 
magnetic QCP of another heavy-fermion system 
Ce(Rui_xRhx)2Si2 = 0, 0.03) using single-crystalline 
neutron scattering llisll . Stoichiometric CeRu2Si2 is an 
archetypal paramagnetic heavy-fermion with enhanced 



C/T 350 mJ/K^ mol and 7k 24 K Ud]. Extensive 
neutron scattering studies of CeRu2Si2 tV/k have shown 
that spin fluctuations possessing three-dimensional 
(J = 3) character are exceUently described by the SCR 
theory for heavy fermions 1^. A smaU amount of Rh 
doping, X > Xc — 0.04 fT§\, induces an antiferromag- 
netic phase (see the inset of Fig. 13 of the sinusoidafly 
modulated structure with the wave vector ^3 = 0.35c* 
LI 9,1 . Samples nearly tuned to the lowest concentration 
QCP (x - Xc) show non-divergent C/r (T 0) |20] and 
Ap oc T^/^ 1T5], which are consistent with the SDW QCP 
in J = 3. Thus one can expect that Ce(Rui_xRhx)2Si2 
(•^ ^ -^c) is suited to investigate the SDW QCP without 
disorder effects. 



EXPERIMENTS AND ANALYSES 

Neutron- scattering measurements were performed on 
the triple-axis spectrometer HER at JAERI. It was oper- 
ated using final energies of = 3. 1 and 2.4 meV provid- 
ing energy resolutions of 0.1 and 0.05 meV (full width 
at half maximum), respectively, at elastic positions. Sin- 
gle crystals with a total weight of 19 g (x = 0) and 17 
g (x = 0.03) were grown by the Czochralski method. 
Two sets of multi-crystal samples aligned together were 
mounted in a He flow cryostat so as to measure a {hOl) 
scattering plane. All the data shown are converted to the 
dynamical susceptibility Imxiq^E). It is scaled to ab- 
solute units by comparison with the intensity of a stan- 
dard vanadium sample. We note that a new point of the 
present work is unprecedented experimental accuracy in 
determining the critical exponent using large samples 
and long counting time. This has enabled us to determine 
the singularity of the QCP and to make qualitative con- 
clusions of the universality class, which was very diffi- 
cult in the pioneering work using the related compound 
Cei_^La;,Ru2Si2 LS,^. 



Previous Study of CeRu2Si2 

Let us first make brief comments on our previous 
neutron scattering study of CeRu2Si2 lHvl] . We showed 
that spin fluctuations of CeRu2Si2 are reasonably well 
described by the SCR theory for heavy fermions 1 6] . This 
result has the following two implications in connection 
with the SDW QCP. First, the spin fluctuations at low 
temperature T = 1.5 K can be parametrized by the SCR 
form 10] 

X{q^E)-'=XL{E)-'-J{q). (1) 

This equation means that the local dynamical suscep- 
tibility Xl{E) = Xl/{^ - i^/rL), expressing the lo- 




FIGURE 1. Observed (a) and calculated (b) intensity maps 
of constant-£^ scans taken with E = \ meV at T = 1 .5 K for the 
sample with x = 1 17]. They are shown on the surface of the 
irreducible Brillouin zone. 



cal quantum fluctuation by the Kondo effect, is modu- 
lated by the intersite exchange interactions [J{q) = 
Y.r^O'^r^o^wi^q ' ^)]- A number of constant- 2 and -E 
scan spectra can be reproduced by Eq. ([l} with the ad- 
justable parameters of Xl, Tl, and 14 exchange param- 
eters 1 17]. In Fig. [I] we show observed and calculated 
intensity maps of constant-^ scans with E = I meV at 
r = 1.5 K. One can see that the calculated intensity re- 
produces the complicated antiferromagnetic spin fluctu- 
ations with the three peaks and weaker structures around 
the Z and N points. 

Second, the temperature dependence of the spin fluctu- 
ations of CeRu2Si2 illv can be approximately described 
by the SCR scenario fa], in which the T dependence of 
Eq. ([U is brought about by the single T dependent pa- 
rameter Xl{T). The temperature dependence of Xl{T) is 
determined by the self-consistent relation |6]. The exis- 
tence of this single T dependent parameter indicates that 
the underlying mechanism is controlled by a neighbor- 
ing SDW QCP. In fact, the numerically calculated Xl{T), 
showing Xl{T) oc const — T^^^ 1 17] in a low temperature 
range, agrees with the computation [Eq. (|2|)] by the renor- 
malization group theory of the SDW QCP |21]. Along 
this line, the purpose of this work is to check whether 
the r^/^ dependence of really occurs in CeRu2Si2 
and in the nearly tuned sample Ce(Ruo.97Rho.o3)2Si2- 



QCPofCe(Rui_;,Rhj2Si2 

In the scenario of the SDW QCP in J = 3, which 
was well established by the renormalization group theory 
lIslBEIl. the wave- vector dependent susceptibility for 
the tuned sample (x = Xc) diverges as x{^3) T~^^^ 



or the characteristic energy of spin fluctuation, i.e., 
the inverse correlation time, depends on temperature as 
r(*3) Xih)'^ T^^^' By taking the detuning effect 
(x < Xc) into account, the leading two terms of Ff^^) 
computed by the renormalization group theory ||4|, |2l|] 
are given by 



^(^3) = Cl + c2^3/^ 



(2) 



where ci and C2 are x dependent constants. This equation 
is an approximation in the temperature range Tfl <C T <c; 
7k, where Tfl is a crossover temperature below which 
the system shows the Fermi liquid behavior fl EHl . 

The imaginary part of the dynamical susceptibility at 
Q = k^^q with small \q\ and is approximated ll2ll] 
by 



\mx{k3 + q,E) 



(3a) 



X{k,)T{h)E 

Dc[Kl + ql + F{ql + ql)]m 

[a quadratic expansion of Eq. ([T}] where Dc and F are 
T independent parameters, and Kc is the inverse corre- 
lation length along the c-axis. In Eq. (l3all the product 
X{k^)T{k^) is theoretically T independent. The two pa- 
rameters Dc and F were determined by using constant- 
Q and -E scans for both samples with x = |17] and 
0.03 at T = 1.5 K. These data were fit to Eqs. ^ con- 
volved with the resolution functions. In Fig.|2|we show 
constant-£^ scans through the antiferromagnetic wave 
vector Q= (101) — ^3, and the fit curves for the sample 
with X = 0.03. The good quality of fitting indicates that 
Eqs. ^ well describe the experimental data at T = 1.5 
K. We obtained x independent values of the parameters 
Z)^ = 98 ± 4 (meV A^) and F = 0. 12 ± 0.01. 

The temperature dependence of F(^3 ) =DcK^ for both 
samples has been determined by performing constant- 2 
scans taken at j2 = ( 101 ) — ^3 . These scan data were fit to 
Eqs. Q convolved with the resolution functions, where 
there are two adjustable parameters F(^3) and xi^s)- 
Several fit results of the constant- 2 scans for x = 0.03 
are shown in Fig.O where one can see that the quality of 
fitting is excellent. We also checked the T independence 
of the parameters Dc and F by comparing the constant- 
E scans in Fig. |2| at T = 4 and 8 K with those calcu- 
lated using the T dependent F(^3) midxiks) determined 
by the constant- 2 scans. The calculated curves in Fig.|2| 
agree reasonably well with the observations. Thus we 
conclude that the theoretical approximation of Eqs. Q 
has been experimentally confirmed, and that the fit pa- 
rameter F(^3) has been determined very precisely. 

The temperature dependence of F(fc3) is shown in 
Fig. |4|by plotting data as a function of T^/^. At low 
temperatures the observed data clearly agree with the 
linear behavior of Eq. Q. In fact, by least squares fit- 
ting we obtained: r(it3) = (0.67 ± 0.01) + (0.0095 ± 
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FIGURE 2. Constant-^ scans taken with E = 0.2, 0.4, 0.6 
and 1.2 meV along the Q = (//, 0,0.65) (a) and (1,0,L) (b) 
lines for the sample with x = 0.03. Data at T = 1.5 and 4 
K are shifted by 0.25 and 0.1 emu/mol Ce, respectively for 
clarity. The curves are calculations using Eqs. (|3j, corrected 
for resolution functions with the same fit parameters as those 
shown in Fig.|3l 
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FIGURE 3. Constant- 2 scans measured at the antiferromag- 
netic wave vector Q = (101) — ^3 for the sample with x = 0.03. 
The curves are fit results using Eq. dial with two adjustable pa- 
rameters, r(^3) and xi^s)- The inset shows the phase diagram 
and Tk of Ce(Rui_;,Rhj2Si2 LiiiSJ. 



0.0021)ri-^^^0-0^ (in units of meV) in the range 1.5 < 
r < 16 K for the sample with x = 0, and r{k3) = 
(0. 129 ± 0.007) + (0.020 ± 0.003)ri-49±0 07 ^^^^^ 
1.5 < T < 8 Kforx = 0.03. Therefore we conclude that 
the observed critical exponent of 3/2 ± 0.1 is in agree- 
ment with the theoretical value of 3 /2 and, consequently. 
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FIGURE 4. Energy width T(k^) of the Lorentzian form 
[Eq. Oal l is plotted as a function of T^/^. The full lines are fits 
to r(A;3 ) = ci + C2 with adjustable parameters ci , C2, and v in 
the low temperature ranges. The dashed line is the calculation 
11711 using the SCR theory 1 6]. 



approximation for the SCR curve in the T range 2.5 < 
T < 13K(4<r^/^<47). 

In connection with the neutron scattering study of 
CeCu5.9Auo.i 111], it was theoretically predicted tT4ll 
that the locally critical QCP is relevant for the two- 
dimensional spin fluctuation |13]. This theory also pre- 
dicted that the SDW QCP is relevant for the three- 
dimensional spin fluctuation, being in accord with the 
present results. Finally we note that, to our knowledge, 
the present work is the first clear experimental verifica- 
tion of the SDW QCP among single-crystalline neutron 
scattering studies on the QCPs of heavy fermions and d- 
electron systems. Assuming that criticalities of QCPs are 
classified into a limited number of universality classes, 
we expect that the SDW QCP remains to be observed in 
other systems. 

In summary, we have demonstrated that the quantum 
critical behavior of the heavy fermion Ce(Rui_xRhx)2Si2 
is controlled by the SDW type QCP in three space di- 
mensions. The inverse correlation time, i.e., energy width 
r(^3) of the dynamical susceptibility, shows a T^/^ de- 
pendence predicted by the renormalization group and 
SCR theories 



that the temperature dependence of spin fluctuation is 
controlled by the SDW QCP in J = 3. The same exponent 
for both X = and 0.03 samples ensures that randomness 
due to Rh doping does not affect the criticality. 



DISCUSSION AND CONCLUSION 

The constant ci in Eq. (|2|l is proportional to the the- 
oretical tuning parameter and, hence, c\ oc Xc — x is nor- 
mally assumed |4]. This assumption is consistent with 
the observed values of c\ and the critical concentra- 
tion Xc = 0.04 ± 0.005. On the other hand, the constant 
C2 is theoretically assumed to be weakly dependent on 
X 10]. However we observed appreciable x dependence 
c^ix = 0.03)/c2(x = 0) ~ 2, which may suggest certain 
unknown perturbations. Despite this problem, we think 
that the critical exponent of 3/2, which is determined 
solely by basic characteristics of the system (the space 
dimension d = 3 and the dynamical exponent z = 2), is 
more important and decisive to conclude the nature of the 
QCP 

An advantage of the present neutron scattering study is 
that Eq. ^ holds in a wider temperature range compared 
to those of indirect measurements of bulk properties. Eor 
example, the leading terms of the specific heat C/T = 
Yq — OcT^^^ were shown to have too narrow T ranges 
to be clearly observed 10, 01 • In Fig.|4] the dashed line 
reproduces the SCR computation of r(^3) for CeRu2Si2 
1 17]. Apart from discrepancy of the coefficient C2, one 
can see that the T^^^ dependence of Eq. (|2|l is a good 
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